dd <- datadist(mad_sheet); options(datadist='dd')
jm <- ols(Josh ~ who + g2 + length + year, data = mad_sheet %>% mutate(who = relevel(factor(who), "Josh")) )
cm <- ols(Charlie ~ who + g2 + length + year, data = mad_sheet %>% mutate(who = relevel(factor(who), "Chuck")) )
tm <- ols(Tyler ~ who + g2 + length + year, data = mad_sheet %>% mutate(who = relevel(factor(who), "Tyler")) )
jmm <- ols(Justin ~ who + g2 + length + year, data = mad_sheet %>% mutate(who = relevel(factor(who), "Justin")) )
trib.cols.rms <- tribble(~term, ~new, ~level, ~label,
'who', 'Album picked by', '1', 'Album picked by',
'g2', 'Genre', '2', 'Genre',
'length', 'Length', '3', 'Length',
'year', 'Year', '4', 'Year',
'Charlie', 'Charlie', '5', 'Charlie',
'Justin', 'Justin', '6', 'Justin',
'Tyler', 'Tyler', '7', 'Tyler',
'Josh', 'Josh', '8', 'Josh')
require(tibble)
prr <- rbind(
rms.sum.table3(summary(jm, year = c(1985, 2005), length = c(40, 60), who = "Josh"),
trib.cols.rms, anova(jm), raw = TRUE) %>% mutate(X = "Josh"),
rms.sum.table3(summary(cm, year = c(1985, 2005), length = c(40, 60), who = "Chuck"),
trib.cols.rms, anova(cm), raw = TRUE) %>% mutate(X = "Chuck"),
rms.sum.table3(summary(tm, year = c(1985, 2005), length = c(40, 60), who = "Tyler"),
trib.cols.rms, anova(tm), raw = TRUE) %>% mutate(X = "Tyler"),
rms.sum.table3(summary(jmm, year = c(1985, 2005), length = c(40, 60), who = "Justin"),
trib.cols.rms, anova(jmm), raw = TRUE) %>% mutate(X = "Justin")) %>%
mutate(p = ifelse(term == "who", highc, axis)) %>%
mutate(f = if_else(term == "who", label, facet))
prr %>%
select(X, label, levels, effect, conf.low, conf.high, any_of(c('d.f.', 'P'))) %>%
mutate(across(where(is.numeric), ~round(.x, 3))) %>%
group_by(X) %>%
jgtt(.)